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Abstract 

The traffic of moiecuiar motors which interact through mutual exclusion is 
studied theoretically for half-open tube-like compartments. These half-open 
tubes mimic the shapes of axons. The mutual exclusion leads to traffic jams 
or density plateaus on the filaments. A phase transition is obtained when the 
motor velocity changes sign. We identify the relevant length scales and charac- 
terize the jamming behavior using both analytical approximations and Monte 
Carlo simulations of lattice models. 



1 Introduction 

Biological cells exhibit complex patterns of intracellular traffic: vesicles shuttle be- 
tween different cellular compartments, others travel from the cell center to the pe- 
riphery or vice versa, but also filaments, RNA, chromosomes, and even viruses are 
permanently on the move within cells [1,2]. An extreme example is the long-ranged 
transport of vesicles, organelles, and proteins along the axons of nerve cells [3], which 
can be up to a meter long. All this traffic is based on the molecular motors kinesins, 
dyneins and myosins which move along cytoskeletal filaments [2,4,5]. These mo- 
tors catalyze a chemical reaction, the hydrolysis of ATP (adenosinetriphosphate), and 
transform the free energy released from this reaction into active movements and me- 
chanical work. In the following, we focus on processive cytoskeletal motors which can 
make many chemo-mechanical steps while staying bound to a filament, in particular 
on the large-scale traffic driven by these motors [6] . 

These molecular motors exhibit movements on several length and time scales which 
range between nanometers and millimeters or centimeters (up to a meter in axons) 
and between microseconds and days, respectively [7]: (I) The single step of a motor 
typically has a size of the order of 10 nm, and is generated through the amplification 
of nanometer-sized conformational changes in the catalytic domain of the motor. The 
chemical cycle takes typically of the order of 10~ 2 s, but the actual physical displace- 
ment is much faster, being almost instantaneous on this time scale. (II) On scales of 
one or a few microns, the motors perform active walks along cytoskeletal filaments 
and move in a directed fashion with a velocity of ~ 1 /im/s. (Ill) On larger length 
and time scales which exceed a few microns or a few seconds, respectively, the mo- 
tors perform peculiar 'motor walks', which consist of alternating sequences of active 
directed walks along filaments and passive non-directed diffusion in the surrounding- 
fluid environment upon unbinding from filaments [6,8]. These motor walks are a con- 
sequence of the fact that molecular motors function in a noisy environment, and that 
the motor-filament binding energy can be overcome by thermal fluctuations. 

* Email: Melanie.Mueller@mpikg.mpg.dc 



1 



Theoretical investigations of molecular motors have largely focused on movements 
on the intermediate scale (II), motivated by single- molecule experiments developed 
during the last 15 years which made the direct observation of the active motor walks 
possible, see e.g., [9-12]. These active walks have been studied using a variety of 
modeling approaches such as various types of ratchet models (reviewed in [13-15]), 
chemical kinetics models and discrete Brownian ratchets [16-19] and Brownian net- 
works [7, 19]. Only little has been done on the dynamics of the single step, i.e., on 
movements in regime (I) , because the duration of a motor cycle still exceeds the time 
scale accessible to Molecular Dynamics simulations. 

The movements over large time and lengths scales, i.e. in regime (III), have recently 
been studied in some detail. Our group has introduced a class of lattice models which 
allows us to study the combination of active walks along filaments and diffusion in a 
rather generic way without making assumptions on the underlying motor mechanisms, 
and, at the same time, to describe specific motors by adapting the parameters to 
the measured transport properties [6]. These models also apply to various types of 
'motor particles' which consist of cargo particles with one or several molecular motors 
attached. 

Using these models, many properties of the large-scale motor walks can be cal- 
culated analytically [20,21]. In addition, these lattice models can be used to study 
motor-motor interactions which are easily incorporated. The most obvious such in- 
teraction is the mutual exclusion or hard-core repulsion which leads to traffic jams 
on crowded filaments [6]. In that case, these models constitute a new class of driven 
diffusive systems or exclusion processes: In contrast to driven lattice models studied 
in the context of non-equilibrium phase transitions, e.g. [22-24], and vehicular traf- 
fic [25], the driving, i.e. the active movement is localized to the filaments and coupled 
to non-driven diffusion through the binding/unbinding dynamics. These new driven 
lattice models have been studied quite extensively both by our group [6, 26-28] and 
by several other groups [29-32] . 

One particularly interesting system is a tube-like cylindrical compartment with 
filaments aligned in parallel to the cylinder axis in an isopolar fashion, which mimics 
the geometry of an axon. Since the lattice models for molecular motors are driven 
systems, the boundary conditions are crucial for the resulting behavior. In our pre- 
vious work, we have studied tube-like compartments with several types of boundary 
conditions: periodic boundary conditions [26], closed boundaries [6,28], and open 
boundaries, which are coupled to motor reservoirs [26]. The latter systems exhibit 
boundary-induced phase transitions. 

In this article, we consider a half-open tube which is coupled to a motor reservoir 
at one end, but has a closed or reflecting wall at the other end. We consider the two 
cases of motors moving towards the open and towards the closed end. This geometry is 
inspired by the geometry of the axon [6], which has a closed end at the axon terminal, 
but, at the other end, is connected to the cell body, where the motors are synthesized 
and which therefore acts as a motor reservoir. 

The paper is organized as follows: In sectional we introduce the model, in sectional 
we present the general properties of these systems and discuss the phase transition that 
is obtained when switching the direction of motor motility. We then characterize the 
motor traffic jams in section 0] where we use an adsorption equilibrium approximation 
to obtain analytical results. In the appendix, the same method is applied to the closed 
tube system studied in [6,28]. We close with a short summary of our results. 
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Figure 1 : Motion of motor particles along a filament within a half-open tube (a) and 
the corresponding lattice model (b). At the left end, the system is connected to a 
motor reservoir with fixed bound and unbound densities Pb,in and p u b,in, while the 
right end is closed. Motor particles bound to the filament perform active movements 
with velocity i>b, while unbound motors perform symmetric diffusion with diffusion 
constant D u b- In the lattice models, these movements are described by the forward 
step probability a on the filament and the symmetric hopping probability 1 /6 to each 
neighbour site for unbound motors. In addition, unbound motors which reach the 
filament bind to it with probability 7r ac i, and bound motors unbind from the filament 
through steps to adjacent non-filament sites which occur with probability e/6 per 
non-filament neighbour site. The motors interact via hard core exclusion. 

2 The model 

To mimic transport in an axon, we study the traffic of motor particles in a half-open 
cylindrical tube with a single filament located along its cylinder axis. The tube has 
length L and radius R. The bound motors move along the filament and the unbound 
motors diffuse freely within the cylinder, see figure Q At the right end, the tube 
is closed, i.e. no motors can leave or enter the system. This mimics the synaptic 
terminal of the axon. At the left end, the system is coupled to a motor reservoir 
which represents the cell body and provides fixed bound and unbound motor densities 
Pban and p u b,in, respectively. 

Each motor can be in two states: bound to the filament, where it moves actively 
to the right with velocity Vh, and unbound, where it performs symmetric diffusion 
with diffusion constant D u ^. The filament is taken to lie along the x-eods, so that the 
system is characterized by the bound motor density Ph[x) and the unbound motor 
density p u h(x). As in [26,28], we neglect the variation of the unbound motor density 
p u b with the transverse coordinates y and z. 

The motors can unbind from or bind to the filament. Since the motors are strongly 
attracted by the filament, the binding rate ttq is taken to be large compared to the 
unbinding rate eo- When motors come close to each other, they may interact. In 
our simple model we include only hard core exclusion which prevents motors from 
occupying the same site. In a mean field treatment, this can be taken into account by 
using exclusion factors of the form (1 — pb) or (I — p u b). 

We are interested in the stationary states of the system. Since the right tube 
end is closed, the net current vanishes in the stationary state. Thus, the directed 
bound current of motors, t>bPb(I — Pb), and the unbound diffusive motor current, 
— 0D u b^:/Oub, must balance to give zero total current: 

d 

«bPb(l - Pb) - 0Alb7^/Oub = 0. (1) 

Here the unbound diffusive current has been integrated over the tube cross section. 
The prefactor <p describes the area available for unbound diffusion. For large radii R, 
4> « irR 2 . Furthermore, in the stationary state, the in- and outgoing currents balance 
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Figure 2: (a) Bound density profiles pb and (b) current densities jb for different 
boundary densities Pb.in as obtained from Monte Carlo simulation (<0>) and from the 
adsorption equilibrium approximation (solid lines). The parameters are vt, = —0.001, 
e = 2/300, 7T = 2/3, = 4. 



on any filament site which leads to 
d 

[«b/0b(l - Pb)] = 7r /9 u b(l - Pb) - £0Pb(l - Pub) (2) 

where the small diffusive part of the bound current has been neglected. The change 
of the bound current drives the system out of adsorption equilibrium for which the 
binding term 7Top u b(l — Pb) and the unbinding term eoPb(l — Pub) would be equal, i.e. 

7!"0Pub(l - Pb) = £0Pb(l - Pub)- (3) 

For simplicity, we assume that the boundary densities at the left boundary, ph 7 m 

and 

Pub, in, fulfill this adsorption equilibrium condition J3J). 

In addition to studying equations Q and (0) by analytical approximations, we 
use Monte Carlo simulations of lattice models as in [6, 26] to obtain the stationary 
density and current profiles. In the simulations, the motors are represented by random 
walkers on a cubic lattice with lattice constant I, see figure Qlb). A line of lattice 
sites represents the filament. Motors at the filament sites perform a biased random 
walk, while motors at non-filament sites perform symmetric random walks. The jump 
probabilities per unit time r are a, (3 and 7 for a forward, a backward and no step at 
filament sites, respectively, and 1/6 for steps to each neighbour site at non-filament 
sites. In the following, we use /3 = for simplicity. Unbinding from the filament 
occurs with probability e/6 to each adjacent non- filament site, and unbound motors 
which reach the filament bind to it with probability 7r a d. Exclusion is implemented by 
rejecting all hopping attempts to sites which are already occupied by another motor. 
The hopping rates of the lattice models are matched to the transport properties of the 
motors via r = £ 2 /(6Z? u b), «b = oi£, 7r = 47r a d/6 and eo = 4e/6, see [6,26] and the 
appendix C in [7]. 



3 Steady states and relevant length scales 

When the motors on the filament move to the right, i.e. for > 0, they jam up at 
the closed right end until the tube far from the left boundary is completely filled with 
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motors. In the bulk of a long tube, one thus has bound and unbound motor densities 
Pb = 1 and p u b = 1. Indeed, these density values are a fixed point of equations 
and J2J)- Within a boundary region at the left end, the densities cross over to the 
boundary values pb.in and p u b,in- 

Likewise, if the motors on the filament move to the left, i.e. for Vb < 0, they are 
driven out of the tube at the left end, so that the bulk is left empty with pb = and 
Pub = 0, which is another fixed point of equations and However, a boundary 
region at the left end survives, fed by ingoing diffusion. Typical motor density profiles 
for this case are shown in figure Eta). They display a "jam" region at the left end, 
separated from the "empty" bulk region by a localized domain wall or shock. The 
system is thus characterized by two length scales: the bulk length scale £ on which 
the densities approach their the bulk values and which is independent of the boundary 
densities, and the jam length L* which describes the width of the jammed boundary 
region and which, of course, strongly depends on the boundary densities. 

On tuning the motor velocity Vb, the system displays two phases, a high and a low 
density phase, which are dominated by the closed right end. Taking the bulk density 
as order parameter, one has a phase transition at Vb = 0. As will be shown both in 
analytical approximation and in Monte Carlo simulation, both characteristic length 
scales £ and diverge with the power law 1/vb as the motor velocity Vb approaches 
zero. 

In an experiment, this limit can easily be realized by reducing the concentration 
[ATP] of the motor fuel ATP, since for small [ATP], v h ~ [ATP]. Changing the sign 
of i>b is much more difficult, one possibility is to use motor particles driven by motors 
of two species with opposite directionality, so that one could invert the movement by 
activating one motor species while deactivating the other. If these motor particles 
switch stochastically between the two directions, one could change the sign of their 
average velocity by influencing the switching rates through regulatory molecules. For 
example, tau proteins which strongly suppress movements towards the synapse in 
axons can induce the retreat of vesicles into the cell body [33] . 

For Vb — 0, i.e. when there is no active motion of the motors, the system is in 
thermodynamic equilibrium, and the adsorption equilibrium as described by equation 
@ is valid for all x along the filament. Thus the bulk densities are pb = Pb,in and 
Pub = Pub. in- This case is the only one where the bulk values are dominated by the 
open left filament end. A phase diagram is shown in figure |3 the control parameters 
are the motor velocity i>b and the left boundary density pb,m- 

To examine the bulk behaviour, equations JQ) and J5J) are linearized around the 
appropriate fixed point. In the case of motors moving to the left (vb < 0) with bulk 
values pb = Pub = 0, the linearization leads to an exponential density profile with 

Pb ~ J\fexp(-x/g) and p ub « ( — - — \ ) Pb (4) 



with the bulk length scale 
£ = 2Ax h 

where 



1 — 4/1 + 4 



Aa; ul 



TO 7T £ 



(Aa; u b) 2 _ eo 0Aib 



Az b 7T V h 



(•5) 



Axb — — and Ai u b = \ — — — (6) 

CO V 7TQ 
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Figure 3: Phase diagram for the half-open tube. For motors moving to the left, i.e. 
fb < 0, the system is essentially empty with bulk density pb = 0, while for motors 
moving to the right, i.e. Vb > 0, it is essentially full with bulk density — I. At the 
phase transition Vb = 0, the system is dominated by the left boundary density and has 
the bulk density pb = pb.in- For Vb < 0, the motors form a traffic jam with length 
at the left end, provided that the boundary density Pb,in is larger than 1/2. Likewise, 
for «b > and pb.in < 1/2, a 'jam of holes', i.e. a region with low density, separates 
the (completely filled) bulk of the system from the left boundary. 



are the average walking distance on the filament before unbinding and the average 
diffusion distance before binding of a motor, respectively. The approximation in equa- 
tion (JjjJ is valid for small motor velocities Vb- The exponential behaviour near the 
bulk value is in agreement with simulations. 

In the case of motors moving to the right, the same reasoning leads to an expo- 
nential approach of the bulk values pb = 1 and p u b = 1 with the length scale 



£ = -2Aa;b,h 




(Al ub ,h) 2 7T0 <^Aib 



Ax 



b.h 



£o Vb 



(7) 



where the last expressions are again valid for small motor velocity Vb, and where 
A£b,h = Vb/TTQ and Ax u b,h = \/ 4>D n b Ao are the average walking distance on the 
filament before unbinding and the average diffusion distance before binding of a hole, 
respectively. This result can be obtained directly from (JjjJ by exploring particle-hole 
symmetry, i.e. by substituting the left by the right boundary condition, pb by 1 — pb, 
Pub by 1 — pub, Vb by —Vb and x by L — x, and by exchanging ttq and eo- 

According to equations iJSJ and Q the bulk length scale £ diverges in the limit of 
small Vb with the power law 1/vb- The same power law is obtained both if Vb = is 
approached from above and from below, but with different amplitudes. 



4 The jam region 

In the following, we will focus on motors moving to the left, i.e. on the case Vb < 0. 
The case of motors moving to the right can easily be deduced as above by invoking 
particle-hole symmetry. Left-moving motors jam up in front of the left boundary, 
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leading to a boundary dominated length scale in the system, the jam length L*. In 
the jam region, the motors move slowly and thus "have time" to equilibrate with the 
unbound motors before moving on. This separation of time scales leads to approximate 
adsorption equilibrium 

7roPub(l - Pb) « £oPb(l - Pub) (8) 

at every filament site x. 



4.1 Density profile and jam length 

With the adsorption equilibrium approximation as given by equation JSJ, one can 
eliminate p u b from the current balance equation Q which leads to 



d_ 

dx 



Pb 



V(l- Kf 
L K 



Pb(l - Pb) 



1 -K 



Pb 



(9) 



where the dimcnsionless constants 



K= ^ and V= Lv * 



7T{) 



<j)D 



ub 



(10) 



have been introduced. The desorption constant K is the ratio of unbinding to binding 
rate, while the dynamic constant V is the ratio of the time i u b = L 2 / ((j>D u b) needed 
to diffuse over the filament length L to the time % = L/vb needed to walk this 
distance. Separation of variables, decomposition into partial fraction and subsequent 
integration of equation (JSJ, using the boundary condition /?b(0) — ph,ixn 

leads to an 

implicit equation for the bound density pb- 



The function 



i 

l-K 



j = ^[g(Pb(x),K) 



K\n(p)--\n{l 



p) 



l-K 2 
K 



■In 



1 - K 



(11) 



(12) 



is a continuous, monotonically increasing, and thus invertible function for < p < 1 
and respects particle-hole symmetry which is equivalent to g(l — p, 1/K) = —g(p, K). 
In the jam region, the profiles obtained from equation (|llf) agree nicely with the Monte 
Carlo profiles, while they differ in the shock and the bulk region, where adsorption 
equilibrium is no longer a good approximation, see figure a). In figure Etb), we also 
show the corresponding profiles of the bound current jb = VbPb(l — Pb) which exhibits 
a maximum close to the domain wall between the jammed and the bulk region. 

We therefore use the position of the current maximum as a definition for the jam 
length L*. In the jam region the current is low due to a too high motor density while 
on the other side of the shock there are too few motors, also leading to a small current. 
As the maximal current is obtained for pb = 1/2, the jam length is given by 



T 



i 

v 



9(\> K ) - g(Pb,in,K) 



(13) 



Note that the jam length L* is positive only if pb,in > 1/2. For smaller boundary 
densities < which implies that the domain wall is located outside the tube 
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Figure 4: Jam length as a function (a) of the boundary density /?t>,in (with fixed 
velocity Vb = —0.001) and (b) of the motor velocity Vb (with fixed pb,in = 0.9). Both 
Monte Carlo (<0>) and analytic results (solid lines) from the jam region approximation 
1113( 1 are shown. The other parameters are as in figure El In (a), we also display 
the averaged bound current Jb as obtained from simulations (O) and from adsorption 
equilibrium approximation ((17(1 with its maximum at pjffn • 



which then does not exhibit a motor traffic jam. The jam lengths determined from 
equation 1(13(1 are shown in figure 0] and agree well with the corresponding simulation 
results. As expected, the jam length increases with the boundary density pb.in- It 
diverges logarithmically w L/(VK) ln(l — Pb,in), as pb,in approaches 1 (of course 
this divergence is truncated by the finite system size) . If the velocity Vb approaches 
zero, the jam length diverges as 1/ub and the jammed region spreads over the whole 
system in agreement with the fact that for Vb = the left boundary dominates the 
whole system, as has been discussed in section above. The exponent —1 of the 
divergence for Vb — ► 0, which has been obtained using mean field and adsorption 
equilibrium approximations, is confirmed by the Monte Carlo simulation results, see 
figure Wlb). 



4.2 The case K = 1 

In the special case K = 1, where binding to the filament and unbinding from the 
filament occur with the same rate, the implicit equation 1(111) for the density profile 
can be inverted and leads to the density profile 



Pb{x) = 



1 ~ Pb,r: 
Pb,in 



exp 



(14) 



The jam length L* in this case is given by 



T 



l 

v 



in 



Pb,in 



(15) 
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4.3 The average bound current 

In order to characterize the overall transport in the system we consider the average 
bound current [6] which is defined by 



Jb 



jb{pb(x)) dx. 



(16) 



Since the bound density is essentially zero in the bulk, the bound current jh(ph) = 
WbPb(l — Pb) has large absolute values only in the boundary region near the left end. 
It is thus appropriate to use the jam region approximation of section T4. II to calculate 
the average bound current Jb, which leads to 



Jb 



K 



V {l-Kf 



l 

l-K 



Pb(J) 



1 

l-K 



Pb,i 



(17) 



For large system size L the right boundary density vanishes, and thus the average 
bound current is given by 



Jb 

Vb 



1 



K 



V {l-K) 



1 - 



1 



1 - (1 - K)pb- 



(18) 



Interestingly, as V ~ Vb, the current Jb in the thermodynamic limit does not depend 
on the motor velocity Vb- This is because the jam length J„ decreases as 1/vb, thus 
cancelling the expected Jb ~ «b behaviour. 

The average bound current from equation 1|17|1 agrees well with the current from the 
Monte Carlo simulation, see figure Ufa) . As a function of the left boundary density 
jOb,in) it displays a maximum absolute value at a density p™?*- This density thus 
optimizes the motor transport in the system. It can be calculated via 







dpb,b 



1 

V 



jb (pb{L)) g'{p h (L)) d ^ L - j h (pb>) g'(PbSn) 



dpb,i 



y [jb (Pb(L)) ~ jb (Pban)] .g'(Pb,in)- 



The last expression follows from 

= g {pb{L))— g (p b ,in) 

OpbSn 

which is obtained by differentiating (|11|) . taken at x = L, with respect to pb,i 
g'(pb,in) > for < Pb,in < 1, a current extremum can occur only for 

jb (Pb.in) = jb (Pb(L)) . 



(19) 

(20) 
. As 

(21) 



Since equality of the densities at the left and right end occurs only for a completely 
full or completely empty tube, the condition for extremal current is 



Pban = 1 - Pb(L). 



(22) 



Using this condition in equation (|ll(l with x = L, one obtains the density that ex- 
tremizes the average current from 



V = g(p^)-g(i- P Z:) 



(23) 



9 



The theoretical calculation agrees well with the simulation results, see figure The 
current extremum occurs when the tube is approximately full, i.e. when the jam length 
L* is only slightly smaller than the tube length L. This is consistent with the assump- 
tion that the average current is mainly supported by the jam region, on which the use 
of the adsorption equilibrium approximation for the calculation of the average current 
is based. 

In the limit of large L, one must thus jam up an infinite tube, for which one needs a 
boundary density pb.in = 1. For this case, equation (|23() can be inverted approximately 
which leads to 

l-K 2 -(1-g) 2 

PbM « 1 - K~^e exp 
i.e. the maximum density 1 is approached exponentially for large L. 

5 Summary and conclusions 

We have studied the molecular motor traffic in half-open tube-like compartments and 
characterized the traffic jams which result from the mutual exclusion of motors. In 
particular, we have obtained analytical results for the jam length and for the param- 
eters which optimize the traffic in these compartments. 

However, although the current within these compartments can be optimized, the 
stationary state of these half-open tubes is always characterized by a density which 
approaches either one or zero far from the open end (with a phase transition separating 
these two cases). Long tubes are therefore either essentially filled or essentially empty. 
In axons which are mimicked by this tube geometry, additional regulatory processes 
are therefore necessary to maintain efficient stationary transport. 

These effects are not restriced to the tube geometry but apply generally to systems 
which are coupled to a particle reservoir at one end and which have two types of 
channels 'perpendicular' to this surface, one with unsymmetric directed and one with 
symmetric diffusive transport. 

A The closed tube 

So far, the motion of motors in a half-open geometry has been considered. A related 
system is one with closed boundaries at both ends which has been studied in [6,28]. 
In this case, motors cannot enter or leave the system, and the total number of motors 
within the tube is fixed: 

N= f [p h (x) + # ub (x)] da. (25) 
Jo 

For this closed system, the same arguments as for the half-open system can be applied, 
however the integral constraint (|25(l makes the analysis more complicated. 

Typical profiles for the closed system and motors moving to the right (v\, > 0) are 
displayed in figureOJa). The profiles show the same two characteristic length scales as 
in the case of the half-open system: the bulk length £ and the jam length L*. The bulk 
length £ which is independent of the boundary densities is again given by equation 



K 



Vb\ 



1 + if 2 D uh q 



(24) 
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Figure 5: (a) Density profiles pb for the closed tube for various particle numbers 
N as obtained from Monte Carlo simulations (0) and from the analytical jam region 
approximation 1(26(1 (solid lines). The parameters are v b = 1/300, eo = 2/300, ttq = 2/3 
and 4> = 4. 

(b) Jam length L* (0) and average bound current Jb (O) for the closed tube as 
functions of the number N of motors within the tube. The data points are from 
simulations, the lines display the corresponding results from the adsorption equilibrium 
approximation as given by equations 1)26(1 and 1(30(1 , respectively. The other parameters 
are as in (a). 



Depending on their direction of motion, the motors jam up at one of the closed 
ends. In this region, the motors move slowly, and the adsorption equilibrium approx- 
imation (jSJ can be applied which leads again to an implicit equation for the bound 
density: 

\ = ^[9( Pb (x),K)-9(p b (0),K)} (26) 

with the function g as defined in equation ((12(1 . This equation corresponds to equation 
(111(1 . however, the left boundary density Pb(0) is now unknown. It can be determined 
from the particle conservation constraint ((25(1 . which for adsorption equilibrium can 
be integrated and leads to: 



V 



N 



L(l 



= F(p h (L),K, cb)-F( Ph (Q),K, 0) 



(27) 



with the function 



F(p,K,4>) = -In 



i 

l-K 



1 



1-(1-A> 



K 



1 



2 1 + 01- {l-K)p 



(28) 



Together with 



V = g(p h (L),K))-g(p h (0),K) 



(29) 



(obtained from equation 1(26(1 by setting x = L) one has two nonlinear equations for the 
boundary densities pb(0) and p\,(L), which are needed in equation l(26|) . The resulting 
density profiles agree well with the Monte Carlo profiles in the jam region, see figure 
Eta). 
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The average bound current and its maximum can be calculated in the same way 
as for the half-open tube, leading to: 



J b 1 K 1 1 



(30) 



«b V{1-KY j^- Pb (L) ^-pbCO) 



with an extremum as function of the particle number N at 



3b(pb(0))=3b(f>b(L)) or p h (0) = l-p h (L). 



(31) 



Thus, from equation l|27|) . the maximum occurs for the particle number 



jymax ^ 



-[F(l- Pb (i))~F( Pb (L))]. 



(32) 



L(l + 0) K 



The results for the average current agree quite well with the simulation results, see 
hgure[5fb). Only near the current extremum the Monte Carlo curve is slightly sharper 
than the curve from the adsorption equilibrium approximation which leads to a dif- 
ferent value for the current extremizing particle number i\T max . 

For a very long tube, the current extremum occurs at particle number 



with h(K, 4>) = -±lnK- a-^a+if 2 ) 2+(3+if)<^ and thug approaches 



for large L. 
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